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ABSTRACT 

We report a mass reconstruction of A1689 using Particle Based Lensing (PBL), a 
new technique for Strong+Weak lensing that allows a variable resolution depending on 
the data density and the signal-to-noise. Using PBL we also calculate the covariance 
matrix for the resulting mass map. The reconstruction of A1689 shows a secondary 
mass peak in the north-east direction confirming previous optical observations. This 
indicates that the central region of the cluster is still undergoing a weak merger. We 
have used this mass map to measure power ratios of the dark matter distribution and 
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compared it to the X-ray distribution. We find that the power in the X-ray distribution 
is lower suggesting a smoother and rounder gas distribution compared to the dark matter 
distribution. We fitted an NFW profile to the profile derived from the mass map and 
we find that the lensing mass within 1 Mpc is 1.5 ± 0.33 x 1O 15 M . This is higher than 
the X-ray mass. 

Subject headings: methods: statistical, analytical, data analysis, gravitational lens- 
ing, galaxies: clusters, A1689 



1. Introduction 

Galaxy clusters are the largest gravitationally bound objects in the universe. They produce 
distinct observational signatures (strong emission in the x-ray and sub-millimeter wavelengths, 
distorted background galaxies) providing information that constrain cosmological structure forma- 
tion, dark matter dynamics, gas-physics and galaxy formation. The number of galaxy clusters as 
a function of mass and redshift, i.e the mass function of clusters will provide strong constraints on 



density f2 m , and evolution of dark energy w(z) (et 


■ 


Francis et al. 


2009 


• 


Bartelmann et al. 


2006; 


Eke et al. 


1996: 


Sheth & Tormen 


1999: 


Frenk et al. 


1990: 


Press & Schechtei 


1974). Cluster 


mass 



measurements are usually plagued by assumptions related to hydrostatic equilibrium and dynam- 
ics of the galaxies. Gravitational lensing measures the mass of the clusters without making any 
assumptions about the dy namical state of the cluster. The prim ary bias in this case is from line- 
of-sight projection effects (jSpinelli et al.ll201ll : iGruen et al.ll201ll ) which is small for very high mass 
clusters ( Bahe et "all bo 111 ). 



Galaxy clusters are still in the process of formation and consequently a large fraction of them 
have secondary structure. The presence of substructure reduces the accuracy of techniques that 
assume the entire mass to be associated with a si ngle peak. Hence it is i mportant to quantify 
substructure from lensing and X-ray measurements. IVentimiglia et alJ (120081 ') have fit an empirical 
dependence of mass on substructure parameters and Jeltema et al. (j2008l ) have shown a strong 
correlation between inaccurate hydrostatic masses and presence of substructure from simulations. 

Gravitational lensing not only constrains the projected mass of the cluster, it gives us informa- 



like the cosmic train wreck A520 ( 


Mahdavi et al. 


2007), the supercluster A901/902 ( 


Hevmans et al. 


2008 


Deb et al. 


2010 


) and many others ( 
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cabe et al. 


2009: 


Kiner h Corless 2007: 


Pedersen & Da 




2007 


Cvpriano et al. 


2001: Williams et al. 


1999; 


Alle 


n 19981). Strong lensing (Zitrin et al. 


2010b 


a) 



lets us resolve the core of the cluster in finer detail providing additional information on substruc- 
ture properties at the core. The combination of weak and strong lensing lets us probe the mass 
distribution both at the core and the periphery of the cluster. In this paper we will reconstruct the 
mass map of Abell 1689 using both weak and strong lensing. 
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Abell 1689 is an X-ray luminous cluster with a velocity dispersion as high as ~ [1600| km/sec (jLemze et al 



20091 ) at a redshift of 0.18 and it is one of the richest clusters observed to date. It is well known for 
its spectacular arcs and many multiple images. There have been several lensing studies of A1689 



from v arious facilities ranging from SUBARU, CFHT, HST and many others (jBroadhurst et al 



2005al ) and different gr oups have found a w ide range of values for the measured concentration 
of the cluster (Table 4. iLimousin et alJ 120071 . hereaft er L07). Strong lensing reconstruction using 
both parametric (L07) and non-parametric methods (|Coe et alJl2010l . hereafter C10) have also es- 
timated a hig her concentration (compared to expectation from ACDM simulations) for an NFW 
profile fitting. iDiego et al.1 (|2005l ) have done a multi-resolution non-parametric mass reconstruction 
for A1689 using strong lensing only for the core of the cluster. 

In this pa per we wil l deve l op a Strong+Weak (S+W) lensing analysis using Particle Based 



Lensing (PBL iDeb et al 



2008 



2010) . The advantage of using PBL is the variable resolution 



that can be obtained in the strong (high resolution) and weak(low resolution) lensing regions. 
Additionally the errors in this technique are well understood. This makes calculation of moments 
from the mass map possible. We apply the S+W lensing with PBL to A1689 and compare the 
mass profile and the substructure parameters from lensing and X-ray measurements. 



2. Lensing Basics 



The lensing notation followed in this paper will be similar to IDeb et al.1 (|2008l . |2010| ) . The 
ultimate goal is to map the two-dimensional dimensionless surface mass density k The lensing fields, 
namely dimensionless surface mass density k, the two components of shear 7 1 ' 2 . The components 
of the deflection field a 1,2 is related to the potential ip by, 
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where ipij represents 



dxidx,, 



and ipi refers to 



dip 



Using PBL the k is related to the potential via a linear matrix operation k given by, 

k = Kip 



(5) 



Similar relations exist for the shear and the deflection field. The redshift dependence of the 
lensing fields : k, 7, a is incorporated by multiplying them with the redshift weight factor given by, 
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z(z) - i^k/EfM ( 6) 

Z{Z) ~ {D ls /D s ){z ^oo) (6) 

The redshift weight factor plays an important role for breaking the mass sheet degeneracy with 
strong lensing data. The weak lensing observables are ellipticities of the background galaxies. In 
the weak lensing regime (k < 1) the ellipticity is given by, 

( E .) =9 . = ^W (7) 
(1 - kZ{z)) 

In the strong lensing regime defined by ((1 — n) 2 — I7I 2 ) < 0, the ellipticity (e % ) is given by, 

_ (1 - KZ{z))fZ{z) 

(e) " 5 - WW () 

The strong lensing observables are the angular positions of the multiply imaged sources. The 
difference of the angular distance between a pair of images is directly related to the deflection angle 
via, 



e A - 6 B = a A Z{z) - a B Z(z) (9) 

where A and B represent the two lensed images from the same source. While combining the 
weak and the strong lensing information we use Equations [8] and [9] simultaneously with [5] to 
map the dimensionless surface mass density k(9,z). 



3. Strong and Weak lensing by Galaxy clusters 

In regions with k > 1 rays of light from distant sources take multiple paths from the source to 
reach the observer producing multiple images. There are certain configurations of lens and source 
galaxies that give rise to giant arcs and multiple images making them very rare. The accuracy of 
measuring strong lensing constraints is very high, however since they are rare and occur only in 
the core of the cluster we need weak lensing to constrain the mass distribution in the outer parts of 
the cluster. Weak lensing causes subtle distortion of background galaxies in regions where k < 1. 

The primary challenge in combining strong and weak lensing data is the difference in scales at 
which the various signals dominate. Since galaxies are intrinsically elliptical and the distortion due 
to lensing is at the level of a few percent , several tens of galaxies have to be averaged to extract a 
weak lensing signal. This limits the resolution of weak lensing mass reconstructions to vary from 1' 
for ground based data to 0.5' for space based data. On the other hand strong lensing occurs within 
the central one arcminute region of the cluster. The detailed information about the structure of 
the cluster core will be smoothed out if the combined mass reconstruction is done with the weak 
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lensing resolution. This situation requires a mass reconstruction technique with variable resolution. 
PBL is designed for this purpose, the resolution of the reconstructed mass is a function of the local 
number density of data and the signal-to- noise of the observables. Thus using this technique we 
get a higher resolution at the core and a lower resolution at the periphery of the cluster. 



A combined likelihood function is given by. 

-i _ _2 

"m 1- Km {i> m } J ^rnn \°n l-re„{^„} 



C oc exp 



(10) 



£ ((q^(W) - cx B m))z{ Zl ) - (e A - e B )Y 



i,pairs 



Here n, m refer to the particle index. The first term of the equation is due to wea k lensing on l y. Th e 



covariance in the data arises because of the smoothing procedure described in iDeb et al.l (|2010l ). 
The second term represents the fit between the multiple images and the deflection field at the 
location of the multiple images scaled by the magnification at the image plane. Maximizing the 
liklihood for anyone of these terms is simple. The weak lensing mass ma p has a resolutio n of 0.5', 



and the strong lensing mass map can have a resolution as low as 10" (jCoe et al.ll2010l ) and the 
positions of the multiple images are fit exactly. The Strong+Weak (S+W) reconstruction has a 
variable resolution, the outskirts have a resolution determined by the smoothing scale of the weak 
lensing data and the core has a higher resolution since we fit to the high signal-to-noise strong 
lensing data. The simultaneous fitting of the strong and weak lensing data is complicated since 
there is some ambiguity in the choice of relative weighting between the strong and the weak lensing 
measurements. We weight the strong lensing constraints with the inverse of the effective error in 
Equation [TTJ The advantage of an S+W reconstruction is that it ensures the mass map at the core 
of the cluster is consistent with the outskirts. 

The contribution to the error of the i th pair is given by a combination of the error in redshift 
and the astrometric error in measuring the positions of the multiple images. The error in measuring 
the positions of the images is og = 0.2" (L07). The error a is given by, 

where A9 is the difference in positions for multiple images, ag are the astrometric errors and a z 
are the errors in redshift. 

The astrometric errors associated with the positions of the multiple images is low. However 
the error in the strong lensing mass reconstruction is dominated by Poisson error. Multiple images 
sample the cluster potential at discrete locations at finite number of points. Thus, even if we fit the 
mass at those finite locations very accurately the overall mass distribution will have higher poisson 
error. The PBL reconstruction also has higher residuals (Figure [3]) for the strong lensing images 
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compared to CIO in order to incorporate the weak lensing constraints which have higher errors due 
to intrinsic ellipticity of source galaxies. 



3.1. Covariance of S+W map 

The minimum of the x 2 = ~ 21og(£) (defined in Equation flUl) gives the best solution for 
the potential ip and hence the mass. In order to calculate the error in the mass distribution 
corresponding to an error in the observables we differentaite the x 2 a t its minimum and get, 



(12) 



Here e is the smoothed ellipticity defined in iDeb et al.l (|2010l ) , and 5 X are the residuals to the fit of 
the strong lensing observables. 



The solution for the potential is obtained by inverting the above equation, 

ip = V w i + V s S x , 

where 



and 



Vw = {y^^fC-^ + w^A^A^ ^(GYC-\ 

\ i i / i 

Vs = (sy?) T c- x & +w s ^2 A l ' T A*\ w s A l ' T - 

\ i i / i 



(13) 
(14) 

(15) 



The covariance in the potential is given by, 



Cov^ = VV> T = V w Cov £ vZ + ^Cov^Vf. 



(16) 



Here we have assumed that the terms V w id x Vj^ and V s 5 x e T V^ go to zero since 



(7a) = 



(17) 



The above relationship holds since the shear has m = 2 symmetry and the deflection field has an 
m=l symmetry. Additionally there is very little overlap between weak and strong lensing data 
since the central region of clusters is masked by the presence of bright cluster member galaxies. This 
ensures that there is no correlation bet ween the weak an d the strong lensing data. The covariance 
in Equation (|16p has been calculated in lDeb et al.l (|2010l ). 



The noise covariance matrix for the reconstructed k map is given by, 



Cov K (C) = KCov^K T , 



(18) 
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Here the reconstructed mass is dependent on the scale £ used to smooth the ellipticities. It is 
important to note that the resolution scale in the strong lensing regime is dictated by the separation 
between the multiple images. 



4. Data 



4.1. Weak Lensing 



The weak gravitational lensing measurements used in this st udy were based on de ep archival 
imaging data of A1689 from the Suprime-cam wide-field imager (jMiyazaki et al.1 12002 ) at Subaru 
telescope. Data representing a total exposure time of 1920s in the V— band and 2640s in the SDSS 
i'— band were retrieved from the SMOKA archive system a nd reduced using the standard SDFRED 
reduction pipeline for Suprime-cam data (|Yagi et al.l 120020 . The limiting magnitude of the Subaru 
data are AB magnitude V=26.5 and i' = 25.9 for a 3a detection within a 2" aperture. 



Following iBroadhurst et al.l (|2005al ) , the V — i' color information from the Suprime-cam pho- 
tometry was used to separate gravitationally lensed background galaxies from foreground and 
cluster galaxies by selecting galaxies that have V — i' color > 0.22 magnitudes redder than the 
red sequence formed by early-type cluster galaxies in a color-magnitude diagram. This catalogue 
was augmented by a partially overlapping set of background galaxies, selected to have an esti- 
m ated photometr i c reds hift larger than z = 0.4, calculated from the CFH12k BRI imaging data 



of iBardeau et al.l (|2005l ). The union of these two background galaxy samples, with an i'— band 
S/N > 7, consists of 9943 objects within a sky area of 33'x23'. Our lensing distortion measure 



ments, based on this g a laxy set, made use o f the IMC AT software package (jKaiser et al 



1995 



Luppino fc Kaiser! 119971 : iHoekstra et al.1 Il998h with an implementation tested against the STEP 
simulations of weak lensing data (jHeymans et al.ll2006l ). 



4.2. Strong Lensing 

The strong lensing data used in this analysis is tabulated in L07. In order to identify all the 
strong lensing images a composite image of HST/ACS data in the four bands F475W, F625W and 
F775W and F850LP are necessary. The accuracy of the modeling of the strong lensing constraints 
depends on the knowledge of the re dshifts of the backgroun d sources. The first strong lensing 



analysis for this cluster was done by IBroadhurst et al.l (|2005bl ) (hereafter B05). In B05 about half 



of the multiply imaged galaxies were identified by eye to create a mass model which was then used to 
detect more multiple images. There was misidentification of a few of the multiple images which were 
corrected in L07 and C10 and additional multiple image systems have been identified. Typically 
one (or two in some cases) of the multiple images have spectroscopically confirmed redshifts, the 
other members of that particular system are obtained by fitting a model and ensuring that they 
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have the same color. 



4.3. X-ray Data 



Clusters of galaxies show very strong X-ray emission. Due to the low efficiency of galaxy 
formation, 90% of the baryons are in the form of intergalactic gas. The deep potential well of 
galaxy clusters traps this gas and heats it to X-ray emitting temperatures. The X-ray temperature 
serves as a proxy for the depth of the potential well and hence the mass of the cluster. The X-ray 
observable is the X-ray surface brightness distribution due to free-free emission. It is proportional 
to the integral of r? e where n e is the electron density. In presence of subs tructure clump s can b e 
visible at varying energ i es an d radii. Details on the X-ray can be found in lMorandi et al.l (|2011al ); 
Riemer-Sorensen et al.l (|2009n . 



We use two Chandra X-ray observations (observation ID 6930 and 7289) from the NASA 
HEASARC archive with a total exposure time of approx. 150 ks. The observations have been carried 
out using ACIS-I CCD imaging spectrometer. We reduced these observations using the CIAO 
software (version 4.1.2) distributed by the Chandra X-ray Observatory Center, by considering the 
gain file provided within CALDB (version 4.1.3) for the data in VFAINT mode. Brightness images 
have been extracted from processed event files in the energy range (0.5 — 5.0 keV), corrected by the 
exposure map to remove the vignetting effects: the two final brightness images has been stacked 
together, in order to increase the signal-to-noise ratio. All the point sources have been masked 
out by both visual inspection and the tool celldetect, which provide candidate point sources. 
The background counts have been estimated by selecting regions of the field not contaminated 
from source emissions. The temperature of the X-ray gas is calculated in radial bins by fitting an 
absorbed optically-thin plasma emission model The deprojected te mperature is used to calculate 
the mass profile using the hydrostatic mass eq uation (ISarazin ll988T ). The resulting mass profile is 
commonly fitted to the universal NFW profile ( Navarro et al.l 19971 ) to constrain dark matter halo 
properties. 



5. Substructure Parameters 



Galaxy clusters that are dynamica lly young still in the process of formation. Simulations 



(Evrard et al 



2008 



Kravtsov et al.ll2006l ) indicate that clusters depart from hydrostatic equilibrium 
at varying degrees suggesting the presence of substructure even in the "relaxed" systems. The 
presence of secondary structure leads to inaccuracies in model dependent mass measurements for 
systems which are very far from equilibrium. Hence it is important to quantify the presence of 
substructures. In this section we will outline techniques for measuring substructure parameters 
from observations. We will calculate two parameters: po wer ratios and ellipticity of the cluster 



halo. Simulations have studied the empirical dependence (jVentimiglia et al.l 120081 ) of mass from 



r 
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X-ray distributions on these parameters. 



5.1. Power Ratios 

In order to quantify the substructure from non-parametric mass maps we use the simplest form 
of parameterization. The surface mass density £ is related to the two-dimensional potential: 

VV = 4vrGS (19) 



We do a multipole expansion of the potential (jBuote &: Tsailll995l ) given by 



tp = ip + i>i —z: (a m cos(m0) + 6 m sin(?n^)) (20) 



mr 

m 



where ipo and ip\ are constants and (r,</>) are conventional polar coordinates. Here a m and b m are 
moments of the surface mass density calculated within a circular aperture. The m th moment in 
the x and the y direction is given by, 



a m (r)= / S(f') (r') m cos(m^)dV, (21) 



r'<r 

b m {r)= f E(r / )(r') m sin(m^)d 2 ^. (22) 

J r'<r 

This technique of multipole expansions is directly related to the dynamical state that results 
from fluctuations in the cluster potential. These moments have the following properties. A cir- 
cularly symmetric mass distribution produces a monopole only term. The dipole term vanishes if 
the coordinate system is set to be at the center of the mass distribution. An ellipse contributes 
to even terms only, thus significant contribution to odd terms would indicate presence of substruc- 
ture. These moments are calculated in a circular aperture. This makes sure that the shape of the 
aperture does not produce any bias. 



The powers are given by, 



P = [aohiOR)] 2 , (23) 
1 

2 m 2 r 2m 



p rn = l^-^{a 2 m + b 2 m ). (24) 



We calculate these powers and calculate their ratio in the form P2/-P0J-P3/ Po and P4/P . These ratios 



arevery sensitive to substructure and describe a wide range of cluster morphologies (jJeltema et al 



20051 ). 
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5.2. Ellipticity 

Ellipticity is another measure of secondary structure. Higher ellipticity of halos indicates the 
presence of sub-peaks in the gas and mass distributio n. There is a strong dependence of ellipticity 
on amplitude of mass fluctuations a% (|Ho et al.l 120061 ) . A higher value of as indicates that cluster 
formation has started earlier and hence the measured ellipticity of clusters in the local universe 
would be lower. Clusters are formed through hierarchical merging of smaller dark matter halos. 
Thus at their infancy they have more infalling matter and are more elliptical. As they virialise they 
become more and more spherical. Thus we expect clusters at higher redshift to be more elliptical 
than low redshift clust ers. This has been confirmed by measu ring higher order moments of the 
X-ray gas distribution (jJeltema et al.ll2005l ; iBuote Sz Tsailll995l ). 



The ellipticities are defined as, 



1 ~ A-/A+ 
1 + A_/A+ 



(25) 



where A±(A+ > A_) are the eigenvalues of the moment of inertia tensor / defined in iDeb et al 
(|2010l ). the position angle of the major axis measured north through east in celestial coordinates 
can be estimated through the following equation: 



6 = arctan 



1 !-IU 



7T 

+ 2 



(26) 



5.3. Measuring Substructure Parameters 



Le nsinq 



In IDeb et al.l (|2010l ) we have discussed calculating moments when the covariance of the mass dis- 
tribution is known. Here we will give a brief recapitulation, since the k is correlated, we will use 
singular value decomposition to define a k' that is uncorrelated. 

K > = ( K ) + U t (k - (k» (27) 

where U is the orthogonal matrix from singular value decomposition of the covariance of the k map. 
The mean density is given by, 

<«> = E -°T K l m (28) 

We use this n' as our surface mass density in equations [22] and inverse weight it by the elements 
of the diagonal matrix from the singular value decomposition. Thus we make sure that the weighting 
is done with the uncorrelated k. The the a m ,b m are weighted by W = w/s, where w is gaussian 
with width given by the radial distance at which the measurements are made. The errors are given 
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by: 

<£W) 

The errors in the power ratio are calculated algebraically from these expressions. In case of mea 
suring ellipticities the weight W = k'w/s and the error calculations are similar to Equation 

X-ray 



Following iBuote Canizaresl (|1994l ): iMorandi et al.l (|2011bl ). we calculate the eccentricity e, cen- 



troid and position angle 9 of the X-ray surface brightness distribution. The parameters ob- 
tained from this method have been computed within circular regions of increasing radius R = 
(200, 400, 600, 800, 1000) kpc centred on the centroid of the brightness distribution. To determine 
these parameters from the }th circular aperture of Pj pixels having rtj counts in pixel i, we calculate 
the moments: 

^m,n = jjr Yl n ^ Xi ~ ^ & " W 71 ^ 2 )' ( 30 ) 

^ iGjth bin 

where Nj = X^ejth bin n *' an( ^ ( s > ^) ^ s * ne centroid given by equation fxi^ = and /xo,i = 0, respec- 
tively. We compute e and within circular apertures with increasing radii. Characterization of the 
errors in this procedure has been performed via Monte-Carlo (MC) randomization of uncertainties 
due to Poisson statistics of the surface brightness: we performed 10 3 MC simulations for each of 
the circular apertures. 

For the calculation of the power ratio, we exploit the method as developed by Buote & Tsai 
(1995, 1996). The X-ray background has been accounted for by calculating the moments for both the 
cluster image and the background image (created via poissonian randomization of the background 
as estimated in regions of the field not contaminated from source emissions) and then subtracting 
the background moments from the observed moments. We compute power ratios within circular 
apertures with increasing radii R. 



6. Results 

We have applied PBL as described in § [3] to the strong and weak lensing data for A1689. 
Figure [1] shows the mass map for 7' x 7' = 1570 kpc x 1570 kpc field of view centered on the cluster. 
The resolution of the mass map in the regions containing strong lensing data, i.e the central 4 square 
arcminutes is ~ 17.5". The rest of the mass map has a resolution of ~ 50". This is necessary since 
the SUBARU data has ~ 12 galaxies per square arcminutes. The higher resolution at the center of 
the map enables us to resolve more secondary structure. The advantage of computing the covariance 
matrix is that we are able to calculate physical parameters (with errors bars) related to the mass 
distribution. The mass map in the left panel of Figure [T] shows a secondary peak in the north-east 
direction to the cluster at a significance of 5a. Compared to previous strong lensing only mass 
reconstruction we are able to probe the mass distribution out to a larger radius. The presence 
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of the secondary structure is revealed on addition of the strong lensing data. The optical image 
shows some second group of galaxies at that location. In order to illustrate this we have provided 
a composite map of the lensing, X-ray and optical image (Figure [2]). The optical data is from 
HST/ACS composite of three wavelengths (F475W, F775W and F850LP). We observe a secondary 
peak at a distance of 180 ±7.5 kpc North-East of the main dark matter peak. 




Fig. 1. — A Strong+Weak lensing mass reconstruction of A1689 using Particle Based Lensing. 
The contours of the field of view represent values of n. The color-scale gives the value of k. In 
this figure north is up and east is toward the left. Two clumps of dark matter are detected, the 
secondary clump is toward the North-East of the main clump. 



In order to quantify these secondary structures, we use the centroid and the ellipticity and 
power ratio as a function of radial distance from the central peak. The centroid of the lensing mass 
distribution is RA = 197.873° ± 0.008, DEC = -1.3437° ± 0.0008 and 30 ± 7 kpc away from the 
central BCG and does not change significantly with radial distance. The X-ray centroid is 25 kpc 
away from the BCG. The centroid for the gas at RA = 197.873°, DEC = -1.341° and the dark 
matter are coincident within errors. The ellipticity of the dark matter distribution is consistently 
higher than the X-ray distribution. It rises sharply within the first 200 kpc to asymptote to a value 
of 0.152 ± 0.016. In the inner 100 kpc (28") the ellipticity is only sensitive to one peak and hence 
it is lower, as the secondary peak gets included in the calculation the ellipticity of the distribution 
rises. The X-ray distrib ution is smoother and rounder. This is co nsistent with the theo retical 

This 



results ([Lee fc Sutdl2003l . see Figure 4 and Figure 5) and simulations (lYoshikawa et al. 



2001 



suggests that the gas pressure is isotropic while the collisionless dark matter has anisotropic velocity 
ellipsoids. The position angle for both the X-ray and the lensing does not change significantly with 



Fig. 2. — Optical, X-ray and lensing composite image of A1689. The optical image is composed of 
ACS images from three wavelengths (F475W, F775W and F850LP). The red contours represent the 
dark matter distribution and the cyan contours represent the gas distribution. The optical image 
has a secondary clump of galaxies coincident with a sub-peak in the dark matter distribution in 
the North-East(NE) direction. The X-ray distribution is smoother, though there is an elongation 
in NE direction. 

radius. At 500 kpc the position angle is 11 ± 5.7° for the S+W mass reconstruction and 13 ± 0.67° 
for the X-ray distribution. 

We also calculated the power ratios shown in Figure [5] for this cluster. The lens mass distri- 
bution has higher power than the X-ray for the second (P2/P0) and third (P3/P0) powers. This 
implies that the dark matter distribution is clumpier than the gas distribution. The fourth power 
ratio (P4/P0) is dominated by noise and values represent an upper limit to the measurements. The 
cumulative power ratios decrease as we approach a radial distance of 1 Mpc from the center since 
the substructure gets smeared out at large radius. The zone of higher power represents the clumpy 
central core and the region of lower power the relaxed cluster at 1 Mpc. 

We have done all measurements from Strong+Weak lensing mass map. The dominating source 
of error in the mass map is the weak lensing uncertainties caused by the intrinsic shape noise in the 
background galaxies and their low number density. The error of the entire map is higher than a 
strong lensing only reconstruction of the central core. The map also has variable resolution implying 
that structure is detected at a variable scale. This can be a possible source for systematic error. If 
there are low mass structures present in the weak lensing regime the 50" resolution will smooth it 
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Fig. 3. — Residuals from the fit to the strong lensing data. Most of the observables are fit within 
1.5" . The mass is simultaneously fit to the strong and weak lensing data, thus not all strong lensing 
constraints are fit tightly. 

away thus biasing the power ratio calculations. The effect of such systematics will be studied from 
simulations in future work. The effect of weak lensing dilution should also be low since we have 
removed object with very high ellipticity. 



6.1. NFW Fitting 



We calculate the cumulative mass derived from the Strong+Weak lensing m ass map by integrat- 
ing; th e mass within a given radial distance Figure [H We fit an projected NFW ([Wright Sz Brainerd 
2000) profile to the projected mass and compare it with results from the X-ray measurements. 
The best fit values to the lensing results are shown by the dashed lines. We find the value of 
r 200 = 1-773 =t 0.32 and the concentration c = 5.56 ± 1.26. The mass enclosed within 1 Mpc is 
(1.5 ± .33) x 10 15 MfTy The pr o jected lensing is systematically higher than the X-ray mass calcu- 
lated in lRiemer-Sorensen et al.l (|2009l ). Also the latest calibrations of Chandra data suggest a 10% 
decrease in the measured X-ray mass. This is increase the discrepancy between X-ray and lensing 
masses. 



6.2. Comparison with previous results 

A1689 has been extensively studied multiple wavelengths. CIO summarize the most recent 
results on concentration in Table 2. They have made the highest resolution mass map of the inner 
2' x 2' region using strong lensing only and obtained a clumpy distribution with zero residuals. 
Compared to the CIO mass map the PBL mass reconstruction has a lower resolution at the core 
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Fig. 4. — Dependence of ellipticity of the dark matter (lensing) and gas (xray) halo on radial 
distance. The ellipticity of the dark matter halo rises sharply in the inner 100 — 200 kpc and 
asymptotes to a higher value. 



primarily to incorporate weak and strong lensing in the same mass map. The advantage of an S+W 
map is that we get a mass distribution out to a large radius. For A1689 we have made a twelve 
times twelve square arcminute field of view. Figure [I] is a zoom into the central seven times seven 
square arcminute region. The larger map has enabled us to calculate substructure parameters out 
to a large radius. The advantage of adding strong lensing to the weak lensing data is the detection 
of the secondary structure identified in the NE direction. This substructure is also detected by 
L07 and is coincident with secondary group of galaxies (Figure [2|) in the optical. Presence of 
substructure in the NE direction has also been confirmed with velocity dispersion analysis of more 
than 500 cluster members suggests the existence of s econdary stru ct ure that coincides with a grou 



of galaxies with north-east to the central region (jCzoskd 120041 ). iRiemer-Sprensen et al.l ((2009) 



have als o detected substructure in the north-east direction. Recent SUZAKU observations in the 
X-rays (jKawaharada et al.l |2010| ) suggests anisotropic gas and temperature distributions in the 
cluster outskirts correlated with the presence of the large scale structure galaxies in the photometric 
redshift slice centered around t he clu ster. Some studies have revealed a discrepancy between lensing 



and X-ray mas ses (IPeng et al 



strong lensing (IMorandi et al 



2009), which have been addressed by triax ial modeling of X-r ay and 



2011bl) and b y non-thermal gas pressure (jMolnar et al.ll2010l ). One 



of the peaks f ound by Leonard et al.l ( 20101 ) is also coincident with the secondary structure in the 



NE direction. Haines et al 



(|2010l ) have observed the A1689 field of view in the K-band and found 
that the light distribution is elongated in the 10 — 12' scale. Valentina Vaccas et al. found that the 
substructure region has a lower radio emission than the overall cluster (jVacca et al.ll201ll ). 
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7. Discussion and Future Work 



In this paper we have reconstructed the mass map of A1689 using Strong and Weak lensing 
data. This is the first non-parametric S+W mass map with a covariance matrix. The mass map 
of A1689 has a higher resolution at the cluster core where we have strong lensing data and it is 
smoother toward the periphery. We have also measured the substructure parameters from lensing 
and X-ray distribution and found that the dark matter is clumpier than the gas distribution. The 
covariance calculation technique developed in this paper applies to grid based \ 2 minimization 
methods for mass reconstruction as well. The substructure analysis can be applied to a sample of 
clusters at varying mass and redshift and test the predictions of ACDM cosmology. Additionally 
we can fit the dependence of s catter in the scaling relat ionships between mass and temperature on 



substructure from simulations IVentimiglia et al.l (|2008l ) and correct the observed scaling relations 
at different redshift. 

The current scenario with observations is ideal for detailed study of cluster properties. There 
are several groups doing multiwavelength analysis of galaxy clusters at different redshift range. The 
Cluster Lensing And Supernova survey with Hubble (CLASH, 0.18 < z < 0.9, 5 x 1O 14 /i~ 1 M < 
M < 3 x 10 15 /i~ 1 M Q .) is observing clusters to produce deep HST images similar to A1689. Future 
surveys like LSST and EUCLID will have lensing data for large parts of the sky spanning a wide 
redshift range will provide a cluster sample with well understood selection function. These future 
datasets will be optimal for the substructure analysis presented in this paper. 
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Fig. 5. — Power ratio as a function of radial distance. The higher power in lensing mass distribution 
than in the X-ray distribution. The power also decreases with distance since the substructure 
dominates at the inner core. 
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Fig. 6. — Comparison of the NFW fit between the Lensing and the X-ray data. The lensing mass 
marginally higher the X-ray mass. 



